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ABSTRACT 


This thesis addresses the development of a novel sea clutter module to enhance the digital 
image synthesizer (DIS) integrated circuit, so that the false-target image generated will be 
more resistant to the electronic protection used within inverse synthetic aperture radar 
(ISAR). The modified DIS design proposed incorporates the DIS false target with a sea 
clutter effect so that the interaction between the reflective surface of the target object and 
the background clutter is emulated. Furthermore, the random Doppler shift of the sea 
clutter derived from the novel DIS design can be used to ease the discontinuity in the 
cross-range domain. The DIS false-target image is then closer to that of a real target. The 
signal processing for the sea clutter module and the DIS and the ISAR compression of the 
modulated return are modeled to evaluate the performance of the modified DIS design. 
Such evaluation is based on investigating the discontinuity of the resultant target image in 
the cross-range domain and the correlation between the image and the background clutter. 
Finally, various radar wavelengths and sea states are used to verify the sea clutter is 
reasonably generated by the module under different environmental conditions. Also, 
different false-target profiles that represent the target in motion are used to confirm the 
integrity of the DIS design. 
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I. INTRODUCTION 


A. ISAR AND DRFM JAMMING 

1. ISAR 

Inverse synthetic aperture radar (ISAR) is a well-known technology for 
constructing high-resolution images of moving targets such as ships and aircraft. These 
high-resolution images, usually in the range and the cross-range (Doppler) domain, 
represent the reflectivity function of the target as a function of range and depend on the 
pulse-to-pulse phase coherence of the reflected returns. The image the ISAR generates is 
the result of coherent processing of the return echoes received at various viewing angles 
as the target rotates within the radar’s field of view [1]. 

In order for the ISAR to obtain a clear target image, the movement of the target 
must be significant with respect to the resolving ability of the ISAR. Normally, the 
motion of the target can be further divided into translational motion, the linear change in 
distance along a direction between the target and the radar, and rotational motion, the 
angular change between the target and the radar. In an ISAR range and cross-range 
(Doppler) image, the range resolution is directly related to the bandwidth of the 
transmitted waveform, while the cross-range (Doppler) is obtained from the Doppler 
frequency gradient generated by the relative rotation of the target and depends on the 
number of pulses integrated. 

The 2-D image consists of resolution cells that contain estimates of the magnitude 
and position of the target scattering surfaces both range and cross-range (Doppler). The 
change in translational motion can be compensated for in the ISAR image processing, 
once it is determined. The Doppler frequency produced by a range-resolved scatterer is 
proportional to the rotation rate of the target and the cross-range distance between the 
scatterer and the center of the target rotation. The orientation of the range-Doppler image 
is determined by the target rotation relative to the ISAR. The range dimension is oriented 
along the radar line-of-sight (LOS) [2,3]. 
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Each reflective return is correlated with a range independent reference function 
that is detennined by the transmitted chirp wavefonn so that the range infonnation 
contained within the pulse is resolved [1]. The range-resolve pulse is then processed 
through an azimuth compression to further extract the Doppler information to construct 
the range-Doppler image. The easiest approach to perform azimuth compression is to 
perform a Fourier transform on the range-resolved pulses. In an actual system, the 
azimuth compression is somewhat more complicated and includes estimating the 
rotational motion during data collection and then calculating the compression coefficients 
that give the sharpest focus. 

ISAR imaging has become more critical in such global military operations such as 
target classification, recognition, and identification. Systems like the U.S. AN/APS- 
137B(V)5 radar and the Russian Sea Dragon maritime patrol radar are capable of 
adopting an ISAR imaging mode to detect, classify, and track targets over sea and air 
[4,5]. They provide the target’s range, bearing, and positional data with both medium- 
and high-resolution images for display and recording. Also, due to the ISARs capability 
for high-resolution imaging, it is extremely beneficial to incorporate the ISAR into the 
launching program of weapon systems to perform target identification and location 
during the pre-launch phase of engagement for higher aim point accuracy and probability 
of kill. 


2. Digital RF Memory Jamming 

Actions taken to confuse or deceive pre-launch weapon designation and targeting 
efforts are known as counter-targeting techniques. Counter-targeting actions include the 
use of low radar cross-section materials on the exterior of the target object, as well as 
stealth and deception devices. The use of low radar cross-section material reduces the 
energy of the pulse return to the radar. Stealthy maneuvering lowers the probability of 
detection, while deception devices tend to confuse the targeting system. The purpose of 
the actions mentioned above is to prevent the enemy’s weapon systems from executing a 
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targeting process and valid lock-on so that a missile will either fail to be launched or will 
miss its target. Unfortunately, these actions are largely ineffective against wideband 
imaging radars like ISAR [6,7]. 

The development of digital false-target image synthesizer (DIS) circuits to 
counter ISAR systems were proposed in the late 90s and have been studied continuously 
ever since. The implementation of the DIS takes advantage of readily available digital 
radio frequency memory (DRFM) technology in which radio frequency and microwave 
signals are sampled and stored at high speed. Together, all these DRFM characteristics 
make it possible for the DIS to change its coefficients and create different structured 
targets adaptively in real-time [4,6]. Consequently, beyond sampling and storing signals, 
DRFM systems also recreate the intercepted radar signals to deceive hostile imaging 
radar systems. 

To jam or create false targets against the ISAR, the DIS circuit first modulates the 
intercepted signals fed from the DRFM (in the form of digitized samples of the chirp 
pulses in amplitude and in phase) to recreate the scattering effect of a desired false target. 
Then, the DIS outputs those modulated signals through a digital-to-analog converter to a 
radio frequency (RF) module, where the signal is frequency up-converted and transmitted 
to the victim ISAR. Once the ISAR collects the modulated pulses generated by the DIS 
and coherently processes them to form an image, the outcome is a two-dimensional 
image that tends to contain the movement of a false target image. Note that a DRFM 
alone cannot perfonn the complex modulations required to synthesize large false targets. 

B. PRINCIPAL CONTRIBUTIONS 

The principal work done within this thesis is described in the block diagram 
shown in Figure 1. The major tasks completed are described as follows. 

The first task completed in this thesis was to perform an IEEE literature search to 
investigate the ISAR and SAR methods that have been proposed to discriminate the DIS 
false targets from the real targets. These discrimination techniques allow the ISAR to 
protect itself from being deceived by structured false targets that are generated by 
DRFMs. This electronic protection (EP) allows the ISAR to reject the false targets. 
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Included in this investigation are the image discontinuity, pulse diversity, and cross-track 
interferometry techniques. The next task is to develop a sea clutter module that works 
with the DIS and determine how well the new DIS design can defeat the image 
discontinuity EP techniques used by these types of imaging radar. 



Figure 1. Block Diagram of Thesis Objectives. 

The second objective achieved is to study the performance of adding the sea 
clutter module to the DIS to be used against the ISAR that is capable of false target 
discrimination. The fundamentals of DRFM operation are also investigated for a better 
understanding of how the DIS takes advantage of the DRFM architecture to achieve 
deceptive jamming, or electronic attack (EA), through the creation of structured false 
targets. 

The technique involves calculating and generating the correct sea clutter and 
represents an important contribution to the DIS design. The DIS signals are pulse-to- 
pulse correlated with background clutter so that the false target is less vulnerable to the 
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ISAR correlated image process. Meanwhile, the false target is added with random 
Doppler shifts so that the false target becomes more realistic and cannot be discriminated 
by the image discontinuity EP techniques. 

In order to generate the sea clutter, the following task is to investigate the 
statistical models utilized to describe the behavior of the sea clutter, including the 
Gaussian model, non-Gaussian model, K-distribution model and Weibull model. Of all 
the models investigated, the K-distribution model is most appropriate since it is able to 
incorporate both the small and large clutter profile structures. For the purpose of 
generating the K-distributed sea clutter, the set of statistical random processes such as the 
Gaussian and gamma random function are used. Upon creation, the clutter is analyzed by 
examining the Doppler amplitude spectrum and amplitude range profile. The pulse-to- 
pulse correlation is also examined. 

Modeling and simulation of the new DIS design is performed for several 
scenarios and the resulting DIS false-target image is compared with that of the previous 
DIS design to investigate how well the perfonnance is improved. In addition, simulation 
of the sea clutter under various sea state conditions is performed and compared to verify 
that the random process adopted is acceptable. Both the coherent properties and Doppler 
profile are also examined. 

C. THESIS OUTLINE 

Chapter II describes the synthetic image processing of ISAR. 

Chapter III explains of how the DIS generates false targets to jam the ISAR. 

Chapter IV describes the three types of EP that can discriminate the false targets 
from the real ones. 

Chapter V presents the new sea clutter module that is incorporated with the DIS. 
The performance of such a module is also evaluated through software simulations in the 
MATLAB program. 

Chapter VI contains conclusions about this thesis and recommendations for future 
study in related fields. 
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II. IMAGING RADAR FUNDAMENTALS 


A. RANGE AND CROSS-RANGE RESOLUTION 

An imaging radar such as SAR or ISAR is able to provide a high-resolution target 
image, in both range and cross-range dimensions, at long distance and in bad weather. In 
SAR, the target is assumed stationary and the radar is in motion, while in ISAR, the 
target motion provides the changes in relative velocity that cause different Doppler shifts 
to occur across the target [8,9]. 

Both types of radar have the same capability to resolve an image in range with the 
range resolution [5] 



where c is the speed of light, and A is the bandwidth of the radar chirp pulse. The 
transmitted chirp pulse can be expressed as [10] 

s(t) = rect(j) exp(j 7r(2 f c t + Kt 2 )) (2.2) 

where K = A/t m and t m is the modulation period. To resolve an image in cross-range 
dimension, however, SAR has a theoretical cross-range resolution 

= f (2-3) 

where D is the horizontal dimension of the radar aperture. From (2.2), it is clear that the 
motion of moving targets (i.e., ships that generate complicated motion such as roll, pitch, 
and yaw over the sea surface) is not taken into consideration. This can result in a serious 
distortion and displacement from the target’s true location in the SAR imaging process. 
Therefore, SAR is mostly used to recognize stationary objects. 

Unlike SAR, ISAR takes advantage of the target’s rotation to obtain the resolution 
in cross-range domain. The cross-range resolution of ISAR is expressed as [5] 
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, _± 

a aISAR ~ ~ 


(2.4) 


where A is the ISAR wavelength, and y/ is the relative angular motion of the target 
during the ISAR observation time, as shown in Figure 2. 


|[ ISAR 


Figure 2. The Relative Angular Motion of the Target. 

From (2.3), it is shown that cross-range resolution depends on the amount of 
angular rotation of the target during the radar observation time when the pulses are 
collected. The relative angular rotation from each part of the moving target can produce 
different ISAR cross-range, Doppler resolutions, especially when there is a tremendous 
rotation in roll, yaw, or pitch. 

The pitch motion of a ship causes the top of the mast to have a larger angular 
motion than the bottom of the mast or the superstructure. These differences cause 
different Doppler shifts in the cross-range resolution. Along with the conventional range 
resolution, a vertical profile of the target along its length dimension is generated, which 
allows the ship to be imaged. Roll motion provides the information about the width of the 
ship. Compared to the length of the ship, the width information is usually less significant 
for recognition. Yaw motion of the ship gives a plan view of the target [11]. 

As the ship moves, its angular rates of motion in pitch, roll, and yaw direction are 
not generally known. Consequently, the ship length and height cannot be determined 
since the ISAR result is only a range-Doppler image. An example of ISAR image is 
shown in Figure 3. 
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Figure 3. a) Photo of the USS Crockett b) ISAR Image of the USS Crockett (From [12]). 

Figure 3a shows a picture of the USS Crockett and Figure 3b shows its equivalent 
ISAR image taken by an AN/APS-137. It can be concluded by a comparison of the two 
images that the airborne ISAR was heading directly to the bow of the ship when the 
image was generated. Notice that the length of the superstructure of the ship in both 
images is about the same. Also, notice that two masts shown in the ISAR image can be 
clearly differentiated from the superstructure. The length of the superstructure, the 
number of masts, and the position of the masts are the key characteristics on which an 
experienced operator relies for ship recognition [11]. 
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B. 


ISAR IMAGING PROCESS 


1. Range Compression 

Pulse compression is a technique adopted to improve the perfonnance of the 
linear frequency modulated (LFM) ISAR waveform. Both stepped frequency and LFM 
can be used. However, LFM is the most popular. The returning pulse at the ISAR 
receiver is compressed to a shorter pulse width r = 1 / A, so that the range resolution is 
improved by a factor of 77 r , which is the so-called pulse compression ratio. The 
received pulse is detected, sampled, digitized, and then compressed as shown in Figure 5. 

Figure 4 is an example that depicts how an ISAR perfonns the pulse compression 
process in range and in cross-range. 



S(m,n) 

Figure 4. ISAR Imaging Compression Process (From [1]). 


The range profde C r (m,n) for a received pulse is obtained from the cross¬ 
correlation of the received signal with the transmitted signal 

C r (m, n) = I(m, n ) 0 S(m, n ) (2.5) 
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where I(m,n ) represents the signal received by ISAR, S(m,n ) represents the coherent 
reference signal derived from the transmitted chirp pulse, and ® represents the cross¬ 
correlation process. According to the correlation theorem, the discrete Fourier transform 
(DFT) of the cross-correlation is equivalent to the product of the Fourier transform, such 
that [5] 

¥{l(m,n) 0 S(m, n)} = I(k,n) • S\k, n) (2.6) 

To speed the process up, a fast correlation process is used. In (2.5), I(k,n ) represents the 
output of a Fast Fourier Transfonn (FFT) of I(m,n ) as [1] 

I(k, ri) = F {/(«,«)} (2.7) 

and S*(k,n ) is the complex conjugate equivalent of S(k,n ) , the output of FFT of 
S(m,ri) as shown below 

S*(M) = (F{SK«)})* (2.8) 

Then, the range profile can be obtained by taking the inverse DFT of (2.5), which 
in turn, takes the DFT of the multiplication of (2.6) and (2.7), and thus (2.4) can be 
rewritten as 

C, (m,n) = F~' [l(k, n) ■S\k , n )} (2.9) 

At this point, the received signal is said to be range-resolved for each pulse and 
ready to be further resolved in cross-range for which the Doppler information can be 
determined. 

2. Cross-Range Compression 

The azimuth compression is accomplished by Fourier transforming the resolved 
range profile obtained from N pulses. The N p profiles can be arranged in a matrix as 
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C r (m,n) 


C rNp (m , n ) for pulse N p 
C rN (m, n) for pulse N p -1 

C r2 ( m,n ) for pulse 2 
C rl (m, n) for pulse 1 


( 2 . 10 ) 


where each row corresponds to a returned pulse and each column is comprised of profile 
elements from the same range bin [5], 

The process of azimuth compression is basically the FFT of the time history, or 
the columns in (2.9), of range profiles, and can be expressed as 

C(m,k) = ¥{Cfm,n)} (2.11) 

In other words, each time-indexed row expression in (2.9) is transformed to its 
corresponding frequency-indexed expression, and essentially yields the ISAR range and 
cross-range image. 

C. CHAPTER SUMMARY 

The concept of ISAR imaging is briefly presented and explained in terms of its 
range and cross-range resolution. Next, the LFM ISAR image compression process 
model proposed in [1] and [5] is introduced as an example of how the received pulse is 
processed to generate an image. The mathematical expression for both the range 
compression and the cross-range compression process are also included in the discussion. 

In next chapter, the EA techniques used against an ISAR system are examined. 
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III. ISAR ELECTRONIC ATTACK 


A. INTRODUCTION 

In modern electronic warfare, the idea of EA is to employ weapon systems or 
devices to prevent the enemy from utilizing the electromagnetic spectrum. Using 
weapons or devices to destroy the enemy’s equipment is usually referred to as a “hard- 
kill,” while adopting jamming devices to deny or degrade the ability of a radar system to 
obtain target information like range, velocity, and angle and tracking is referred to as 
soft-kill techniques [13]. 

There have been many efforts devoted to researching different types of jamming 
techniques to defeat ISAR high-resolution imaging capability. Most of the methods 
proposed are divided into two categories: noise jamming and deception jamming. The 
major difference between noise jamming and deception jamming is on the correlation 
between the jamming signal and the real target return. 

The objective of noise jamming is to inject an interference signal into the victim 
radar receiver such that the actual return is completely overwhelmed by the interference. 
In order to perform noise jamming, a jammer has to continuously generate noise that can 
be either modulated in amplitude or in frequency [14]. 

At the victim radar receiver, the ISAR return signal is weaker in energy compared 
to the noise. Although the noise signal does not compress like the target signal does, it 
can still distort the image classification process. The loss of information on the target 
affects the ISAR motion compensation process and can also distort the resulting range 
and cross-range image quality. Since random noise is completely decorrelated with 
respect to the compression process within the ISAR receiving system. Noise jamming is 
not as effective as one might expect against ISAR [13]. 

The goal of deception jamming is to mask the real target by sending suitable 
modified replicas of the real return to the victim radar so that those replicas can go 
through the radar’s signal processing and create a false target to confuse the radar. The 
deception signal transmitted by a jammer basically contains certain inter-pulse 
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information of the ISAR returning signal and has more correlation with the return signal 

[15] . As a result, the deception signal can get through the ISAR imaging process and 
generate a realistic false-target image, which leads to ineffective decision-making by the 
radar, radar operator, or even commanding officer. Deception jamming is generally 
considered as an economic and effective method against ISAR compared to noise 
jamming and will be the technology discussed in this chapter and throughout this thesis. 

B. DIGITAL IMAGE SYNTHESIZER 

1. Overview 

The digital image synthesizer (DIS), proposed in [1] and fabricated and tested in 

[16] , is a pipelined, all-digital image synthesizer designed to deceive LFM ISARs of 
various carrier frequencies. The essential concept of operation is to emulate false-target 
signals and transmit those false signals back to the ISAR so that a false target is generated 
on the resulting range-Doppler image at the end of imaging process. 

As an application of deceptive electronic attack, the power of the false-target 
signal should be high enough to cover the true ISAR return and low enough to avoid the 
saturation of the ISAR receiver, which otherwise could possibly reveal the presence of 
jamming. Also, in order to effectively deceive an ISAR, the false-target signals must be 
appropriately generated such that they are realistic in size and in Doppler effects to reflect 
a feasible target size and motion. 

With the help of modern digital signal process techniques, such as high- 
performance DRFMs, it is possible to enhance the false-target signal with proper phase 
and gain modulations by integrating information including range, Doppler and radar 
cross-section (RCS) into the intercepted ISAR pulses before retransmitting. A block 
diagram of the DIS architecture is shown in Figure 5. 
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Figure 5. DIS Block Diagram (From [1]). 

The following description explains the signal processing of the DIS architecture 
throughout this chapter and is theoretically and mathematically based on the discussion 
contained in references [1] and [16]. 

2. Digital RF Memory Phase Samples 

Most of the advanced military radars process signals coherently to preserve the 
phase of the detected signal against the deceptive repeater jammer. The old-fashioned 
non-coherent jamming device is less effective and the need to replace these with newer 
coherent jammer techniques increases. DRFM-based jammers that store and replicate the 
intercepted radar signal coherently are one of the recent methods used to break through 
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the coherent processing of radar and maintain a reliable jamming effect. Figure 6 
represents the block diagram explaining the flow of the DRFM signal processing, which 
is described in the following paragraphs [9], 



Figure 6. Functional Diagram of DRFM. 

The typical DRFM process includes the following steps: (1) down-converting the 
frequency of the intercepted signal to intermediate frequency (IF) with a reference local 
oscillator (LO) and decomposing it into the in-phase (I) and quadrature (Q) components, 
(2) sampling and digitizing using an analog-to-digital (ADC) converter, (3) storing the 
sampled pulse information into the memory for generating the jamming signal purpose, 
(4) reconstituting the I and Q components modified with the jamming signal to analog 
form through a digital-to-analog (DAC) converter, and (5) up-converting the jamming 
signal with the original LO [6,8,13]. 

Among the processes mentioned above, the conversion of signals into a digital bit 
stream is limited by the speed of the ADC; as a result, it is necessary to choose an IF that 
matches the converter sampling rate. Also, the two oscillators must be very stable so that 
the phase information of the signal is preserved. The signal samples stored in the memory 
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can be utilized for various electronic counter-measure applications, including the 
generation of noise, synthesis of a different in-band frequency, and modulation of false 
targets. 

For example, considering a LFM ISAR that transmits an up-chirp pulse, the 
complex envelope of the transmitted signal can be represented as 

s(t ) = rect(—)e jnAt /r (3.1) 

T 

and that of the pulse intercepted by DRFM is 

Si =rect(-)e J2 ^ fdt+A,2l2T) (3.2) 

r 

where rect(t / t) is the rectangular function, r is the uncompressed pulse width in 

seconds, A is the modulation bandwidth in Hertz, f d is the Doppler frequency due to the 

closing rate between the ISAR and the platform with the DRFM. Notice that, in (3.2), f d 

is usually in tens of Hertz and can be ignored if compared to A, which is in megahertz 
(MHz). Thus, it is reasonable to assume that the phase change within the chirp pulse is 
constant and (3.2) can be replaced with (3.1). 

Depending on the type of DRFM used, the signal expressed in (3.1) can be phase- 
sampled or amplitude-sampled at a sampling time given by 


m 


t =-1 -n- PRI 

fs 


(3.3) 


where m is the index of the intrapulse samples, n is the index of pulse number, f s is the 

ADC sampling frequency, and PRI is the ISAR pulse repetition interval (PRI) in seconds. 
Although many imaging radars use variable PRI, the PRI must stay constant within the 
time to construct one image. Substituting (3.3) for the sampling time,?, in (3.1), the 
sampled phase of the DRFM intercepted signal can be expressed as 

</> 0 (m,n) = —(^ + n-PRI) 2 (3.4) 

r f s 
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In the case of a ^ -bit phase-sampling DRFM, the phase information in (3.4) is 

generated with k p -bit resolution by 2 kp integers ranging from 0 to 2 kp -1. Thus, after the 

DRFM sampling processing, the phase information expressed in (3.4) is stored in the 
memory in the following fonn 


n ) 


M m >n) 

2n / 2 kp 


(3.5) 


where |_xj represents the greatest integer smaller than or equal to x . The resulting phase 

information is then inserted into the DIS process. The identical concept of quantization is 
also applied to a k q -bit amplitude-sampling the DRFM. However, it is necessary to 

convert the ^ -bit amplitude samples into phase before they are inserted into the DIS 
process as shown in Figure 6 [5,16], 


3. False-Target Image Profile 

In order to synthesize the false-target image with the phase and amplitude 
information input from DRFM (both the I and Q channels), appropriate modulations must 
be performed at each range cell to emulate the Doppler shift due to target rotation and the 
expected radar cross-section value of the reflective surface at each location in the range- 
Doppler plane. It is the microprocessor block in Figure 5 where a false-target image 
profile and the corresponding phase and gain coefficients used to synthesize the false- 
target image are generated. 

A false-target profile is constructed by creating an TV,.by N d matrix of scattering 

elements, where TV,, is the number of the resolution cells in range and N d is the number 

of Doppler resolution cells in cross-range. An example of false-target profile for a ship is 
shown in Figure 7. 
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Range resolution cells (r) 


Figure 7. Range-Doppler Profile Template for A False Target. 


For each ( r,d ) cell in Figure 7, its amplitude, A(r,d), depends on the radar RCS 
of the scatterer contained within the cell and according to the radar equation, this RCS 
value is related to the distance between the ISAR and the scatterer, R , in the following 
form 

i 

R qc RCS 4 (3.6) 

The phase of the scatterer within each (r,d) cell, on the other hand, depends on 
the Doppler shift, /(r, d ), due to the relative angular rotation of the scatterer along the 
LOS of the false-target motion. Both A{r,d) and f(r,d) are varied on a pulse-to-pulse 
basis. Therefore, each cell in the profile can be expressed in a tri-dimensional matrix 
comprised of range cell, Doppler cell, and the number of pulses, as follows 

T(r,d,n) = A(r,d)e- J2 * nr ’ d)nPRI (3.7) 

where n is the index of pulse number. Furthermore, the signal of the false-target profile 
for each pulse at each range cell, r , can be obtained by summing all Doppler components 
within each r and is expressed as follows 

Nd(r) 

T'(r,n)='Y j T{r,d,n) (3.8) 

d =1 

The phase of each false-target return is calculated as 
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(3.9) 


tf> T (r,n) = angle {real \T (r,n)} ,imag\T (>,//) jj 

where the angle function is equivalent to taking arctan function on the division of the 
imaginary and real components of T ( r,n) . Once the false-target return and its 
corresponding phase are determined, the DIS phase and gain coefficients, </> inc {r,ri) 
and g(r,n) , are independently generated for each range cell. 

The phase coefficient must be properly generated so that it is recognized by the 
ISAR as a successive phase variation on a pulse-by-pulse basis. Hence, for pulses after 
the first one (n ^ 0), the modified phase coefficient for the range cell r , <f> inc {r,ri) , is 

obtained by subtracting the phase calculated for the current pulse from the sum of the 
phase and the phase variation calculated for the previous pulse. The mathematical 
expression is written as 

<tine (a n) = </>' mc (r, n- 1) + <j) T (r, n- 1) - 4 (r, n) (3.10) 

As to the first pulse ( n = 0 ), it is assumed to remain in phase with no variation such that 
the modified phase coefficient for the first pulse is 

4(h0) = 4(f0) (3.11) 


The phase coefficient in (3.11) also needs to be quantized to an integer value in 
accordance with the k p -bit phase-sampling DRFM. Thus, the modified phase coefficient 

is given by 


4c(a») 


(t>' in c(r,n)2 k " 

2n 


(3.12) 


The modulation of the amplitude coefficient is similar to that of the phase 
coefficient. There are various methods that can be used to produce the gain coefficient, 
g(r,n ), such as unifonn, exponential, and log thresholds. The exponential scheme is 
used for the DIS design. In Figure 7, the gain block multiplies the complex signal by 
2 g(r ' n) , where g(r,n ) represents the number of left-shifts applied [1,16]. Such gain 
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modulation is done by shifting bit locations; the use of multiplexers is eliminated for 
better performance. A detailed discussion on the implementation of the multiplexer is 
found in [16]. The first step to obtain the gain coefficient is to normalize the magnitude 
of T ( r,n ) in (3.8) so that all magnitudes have a range of values between 0 and 1 

T N (r,n)= 1 ■ ' ■ (3.13) 

max i (r,n)\ 

The normalized magnitude is then quantized into 2 />s levels, where k g represents 

the number of bits controlling the left shifts. Due to the hardware limitation, some bit 
values may produce the same amount of shift, which results in fewer control codes 
available than 2 kg . For instance, a gain shifter with a 4-bit control produces only 11 
unique shifts instead of 16. Table 1 shows the possible quantization for the gain shifter 
implemented with a 4-bit control. 


Normalized Magnitude 

Gain Coeffcient g(r, ri) 

Effective Gain 2 s(r '" ) 

0.8-1.0 

10 

1024 

0.4-0.8 

9 

512 

0.2-0.4 

8 

256 

0.1-0.2 

7 

126 

0.05-0.1 

6 

64 

0.025-0.05 

5 

32 

0.0125-0.025 

4 

16 

0.00625-0.0125 

3 

8 

0.0032-0.00625 

2 

4 

0.0016-0.00625 

1 

2 

0-0.0016 

0 

1 


Table 1. Gain Modulation Quantization Table (From [5]). 
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4. The DIS Signal Processing 

As shown in Figure 6, the DIS processes the k p -bit phase samples in both the I 

and Q channels from a phase-sampling DRFM and uses a high-perfonnance 
microprocessor for false-target image control. Each of the DRFM sampled phases 
</>(m,n ), as shown in (3.5), is entered in parallel to N r range bin modulators that are 

derived from the N r range resolution cells of the false-target profile generated by the 
microprocessor. Along with the phase and gain coefficients, every <f>(m,n) is processed 
through a phase adder, a look-up table (LUT), a gain block, and a summation adder 
within each range bin modulator. 

From Figure 5, the DRFM sampled phase is first added to the phase coefficient 
</> inc (r,n) through the phase adder and the resultant output is given as 

$(r,m,n) = 0(m,n) + 4 inc (r,n) (3.14) 

which implies the phase of the m‘ h sample of the n th pulse is added with the phase 
coefficient at the r th range cell of the same pulse. 

The k - bit phase adder output <fr(r,m,n) is used by the LUT to index a sine and 

cosine value to construct a complex signal of unit amplitude with I and Q components, 
respectively. Thus, the LUT output is 

L(r,m,n) = cos(0(r,m,n)) + j sin(0(r,m,n)) (3.15) 

where each of the I (cosine) and Q (sine) output channels has amplitude values ranging 
from -1 to 1. 

The LUT output (both I and Q channels) are then amplitude modulated at the gain 
block in each range cell modulator. The effective gain derived from the gain coefficient 
on which a power of two’s function is applied is referred to as the modified RCS for the 
false-target image that is going to be returned to ISAR. The gain modulated channels can 
be expressed as 
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S(r,m,n) = 2 s(r ' n) L(r,m,n) = 2 g{r ’ n) e mm ’ n)+ ^ r ’ n)) 


(3.16) 


Following the gain modulation, all further range cell processors are summed with 
the current processor and passed forward. The final adder block within each range cell 
modulator sums the result after amplitude modulation with its adjacent (delayed) adder 
output and sends the result forward to the adder in the next modulator. The final complex 
output pulse is collected at the adder in the first range cell modulator. In other words, 
each output pulse is the superposition of copies of the pulse, each delayed with respect to 
one another by the adder delay, scaled differently by the gains and phase rotated by as 


I(m ri) = 


+ z~ 


+ z~ 


+ z 


[ 2 g(1 " 

T 2 

I 2 


n) e j(0(m,n)+^ m (\ 


2 r^g(2,«) „j(<Km,n)+f nc (2 


»)) J 

H)) ]+... 


(N r - 1) f ^g(N r -l,n) j(<f>(m,n)+<f> inc (N r 


-!,«)) J 


(3.17) 


where z 1 represents the delay inherent in each adder[l]. 

For simplicity, the final output of the DIS false-target image signal for the m th 
sample of the n‘ h pulse is rewritten as 


I(m,n ) = S(r,m,n) 2 8(r,n) e 


j(0(m-r,n)+<f> inc (r,n)) 


(3.18) 


r=0 


r =0 


Finally, the digitized DIS output signal is passed to a DAC and transformed into an 
analog signal. The same LO used during the ISAR signal reception phase is further 
applied to up-convert the output signals onto the ISAR carrier frequency. The final up- 
converted jamming signal is then retransmitted back to the ISAR for image processing, 
which is discussed in Chapter II. 


5 Variant of the DIS Design with Multiple False Targets 

everal efforts have been devoted to improving the DIS design since it was first 
proposed. In addition to improving the software codes, algorithms for the reduction in 
response time, and a better perfonnance in parallel processing operations, other efforts 

have further modified the phase coefficient so that the original DIS design is able to 
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generate multiple false targets. In [17], an algorithm that generates multiple false targets 
is proposed, and an example of the resulting ISAR image is shown in Figure 8. 



Figure 8. ISAR Image with Multiple False Targets (From [17]). 

In order to accommodate the original DIS design to generate the effect of multiple 
false targets on the ISAR image, it is the phase of the false-target return at each range cell, 
(/> T (r,n ), in (3.9) that must be modified as [17] 

</>T_ mo <i{ r ’ n ) = kJ T {r,n) (3.19) 

where (j) T mod is the modified phase and k m is a phase modulator with the following form 

*,.=(-1)' (3-20) 

The phase increment in (3.10) becomes 

kc (a n) = (f inc (r, n - 1) + <f ) T _ mod (r, n -1) - <f> T mod (r, n ) (3.21) 

From (3.20), it is clear that the phase modulator reverses the phase by 180 degrees at each 
of the odd range cells. In other words, half of the range cells in the false-target profile 
have opposite phase values while those of the other half remain unchanged. 
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Compared to (3.17), the final false-target signal becomes [17] 


I(m, n) = 1 2 ^( r ’ W ) gj(t( m > n ) + 4inc ( r >«)) 

r =1 

N r /2 

_ gj</(m,n) I 2 g(2r,n) gj<f> inc (2r,n) _| _ gjf(m,n) ^ 1 


AC/2 


r =1 


r=l 


(3.22) 


where two identical false targets are generated. The image containing the term e j<l>inc(2l ' n) is 
attributed from the even-range cells of the false-target profile and the other with the term 
e -M«( 2r - 1 .") j s f rom the odd-range cells. 

C. CHAPTER SUMMERY 

A discussion of the electronic attack measures against ISAR, both the noise and 
deception jamming, is provided in the beginning of this chapter. Next, the DIS concept is 
mathematically introduced, including the signal processing at the DRFM, the generation 
of false-target profiles and the final false-target signal. Finally, a multi-false-targets 
jamming technique, modified from the DIS design, is presented, and the difference 
between the original and new DIS algorithms is described mathematically. 

In the next chapter, the EP techniques that can be used by the ISAR to protect 
itself from deception EA are described. 
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IV. METHODS TO DISCRIMINATE FALSE TARGETS 


As mentioned in Chapter III, jammers that use DRFM-based DIS false targets are 
considered an effective scheme to confuse SAR or ISAR for target image acquisition. As 
a result, the capability of recognizing whether an imaged target is real becomes more and 
more critical for radar reliability. 

There are several methods or algorithms developed to discriminate false targets 
from real targets and to filter them out from the operator’s view. Methods that have been 
investigated include those based on the target’s continuity, pulse diversity, and cross¬ 
track interferometry which represent the three general schemes that have been described. 

A. CONTINUITY OF THE FALSE TARGET IN THE ISAR IMAGE 

The recognition of a DIS false target can be perfonned by inspecting the 
continuity of the target within the ISAR image. An example is shown in Figure 9. Figure 
9a is the image of ship taken by the AN/APS-137 ISAR radar, while Figure 9b is the one 
obtained from the simulation of the DIS false-target based on the size and shape of the 
same ship. 



Figure 9. Comparison of ISAR Images on Real Target and DIS False Target. 


In both images, the superstructure of the ship can be identified and the size of the 
ship can be approximately estimated but it is not difficult to distinguish the real target 
from the decoy. Due to the complex reflectivity of the target’s surface and the Doppler 


27 





shift caused by the background clutter, the ISAR imaging processing cannot perfectly 
compensate for the target’s motion in range and cross-range dimensions [4]. Hence, the 
ISAR image of a real target in Figure 9a is blurred, especially at the edge of the surface. 
Also, the image of the ship’s superstructure is continuous in both range and cross-range 
domains, which means there is a correlation in Doppler shift between consecutive 
reflective surfaces. 

The simulation of the ISAR image for the DIS false-target in Figure 9b, on the 
other hand, is discontinuous in both dimensions. Such discontinuity of the DIS false- 
target image is fundamentally derived from the design of the false-target profile in the 
DIS processing. 

For the DIS false-target profile in Figure 8, the Doppler shift assigned for each 
surface of the false target within a range cell is identical on a pulse-to-pulse basis. In 
other words, for a range cell along the radar LOS, every reflective surface is assumed to 
have the same Doppler shift on each pulse. This assumption implies that there is no 
correlation between consecutive reflective surfaces and results in a gap in the cross-range 
domain among all range cells. The discontinuity in Doppler occurs during the ISAR 
cross-range compression process. 

Another reason the DIS false target is discontinuous is the digitization of the 
phase information in order to meet the DRFM architecture. This digitizing process 
actually transforms the phase information into an integer representation, which causes a 
gap, or an error, between two consecutive cross-range cells at each range cell. Therefore, 
the discontinuity in cross-range is produced after the ISAR cross-range compression 
process is performed [4]. 

Other than taking advantage of the deficiency of the DIS design to determine the 
continuity of the image, altering the ISAR’s Doppler resolution can also be an option that 
enlarges the discontinuity of the target image and helps in identifying the false target. 
Figure 10 shows simulation results of an ISAR image in three different radar Doppler 
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resolutions. These simulations are based on the same DIS false-target profile and can be 
used as an example to illustrate the relationship between radar Doppler resolution and the 
discontinuity of false-target images. 



(a)Doppler cell = 64 (b) Doppler_cell=128 (c) Doppler_cell=256 

Figure 10. Simulations on ISAR False Target with Different Doppler Resolutions. 


From Figure 10a, 10b, and 10c, it is obvious that the finer the Doppler resolution, 
or the more the Doppler resolution cells, the easier to pick up the gap between two 
consecutive cross-range cells. Thus, an experienced radar operator can examine the 
continuity of a target image and determine the authenticity of the target by changing the 
radar Doppler resolution. 

To sum up, the method for ISAR to discriminate DIS false targets by examining 
the continuity of the false-target image is probably one of the most economical methods 
that focus on the flaw of the DIS design. The DIS design requires a modification to defeat 
such EP operations. 

B. PULSE DIVERSITY WITH A DRFM PENALIZING ALGORITHM 

The strategy that discriminates the false target using pulse diversity was originally 
proposed to enhance the performance of SAR against the DRFM-based jammer. It is 
possible to deploy this technique on ISAR to neutralize the DRFM-based jamming signal 
that contains the DIS false-target signal; though few efforts have been devoted in this 
area. Since the pulse diversity method is developed for SAR, the following description 
will be presented from the SAR point of view. 
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Aperture radars utilizing the pulse diversity method are involved in an alternation 
of the transmitted signal and/or its parameters in the fast-time domain at each PRI [18]. In 
most of the cases, a random phase is added onto each of the radar pulses initially, on a 
pulse-to-pulse basis before transmission. Such initial phase information is kept within the 
radar where the DRFM-base jammer has no access. When the DRFM jamming signal 
based on the previous pulse arrives at the radar receiver, its phase information does not 
match with that of the current radar pulse and is mostly likely to be filtered by a matched 
filter. As a result, the radar becomes more resistant against a DRFM-based jammer. A 
DRFM penalizing algorithm is further proposed in [19] and [20] to enhance the pulse 
diversity method. The DRFM penalizing algorithm is a two-dimensional, range-Doppler 
signal processing that effectively negates the DRFM-based jammer image. A system 
block diagram illustrating the process for both the pulse diversity and the DRFM 
penalization is shown in Figure 11. 


Transmitter 
Matched Filter 



Figure 11. Block Diagram for Radar’s DRFM Penalization in Range Process. 


In Figure 11, the signal received by the radar is first processed in sequence 

through a set of matched filters, including the transmitter matched filter and jammer 

matched filter. In order to construct the jammer matched filter, some characteristics of the 

jammer, like its transmitting pulse sequence, must be known beforehand. 
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Consider a SAR system with the following pulse sequence 

P,M m = 1 ,... (4.1) 

where m is the index identifying the discrete measurement in the radar slow-time domain, 
u , and t represents the time index in the radar fast-time domain. The corresponding 
jammer pulse sequence is assumed to be 

m = 1,... (4.2) 

Thus, the signal received by the radar at m‘ h aperture position can be expressed as [20] 

S ( t, u m ) = S ACT (t, u m ) + Sj(t,u m ) (4.3) 

where S ACT (t,u m ) represents the signal returned from the actual target, and Sj(t,u m ) 
represents the DRFM false-target signal. Upon reception, S(t,u m ) is down-converted to 
base band and applied to the transmitter and jammer matched filter, respectively, where 
^(^O and S 2 (t,u m ) are outputs generated. 

The received signal is match filtered with the m th pulse to reduce the effect of the 
jammer signal, and S x (t,u m ) is generated as follows 

S x (t,u m ) = S(t,u m )*P*(-t) 

= S ACT (T U m ) * P* (-t ) + Sj (t, u m ) * P* m (-0 

Since Sj(t,u m ) is the signal based on the information of the previous pulse, which differs 

with the current pulse in phase, the jammer effect within S x (t,u m ) is weakened. An 

example picture of the resultant S x (t,u m ) from [19] is shown in Figure 12, in which the 

strength of the jammer signal is found to be approximately 20 dB lower than that of the 
real target return. Such result can be further considered as the effect of deploying pulse 
diversity in phase against the DRFM false-target jamming. 
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7 7.5 8 8.5 9 9.5 10 


Range, km 

Figure 12. Radar Matched Filtered Signal 5j(7, u m ) (From [19]). 

The received signal is also applied to the jammer matched filter that is based on 
the previous pulse and generates an output, S 2 (t, u m ), with the following expression 

S 2 {t,u m ) = S(t,u m ) * P* ,(-0 5) 

= S ACT & U ,n ) * K-l ("0 + S J (P U ,n ) * P L ("0 

By the same token, the real target return signal, S ACT (t,u m ), is weakened due to the 
mismatched phase between the current and previous pulses. An example picture of the 
resultant S 2 (t,u m ) from [19] is shown in Figure 13. 
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Range Reconstruction Using Jammer Matched Filter 



Figure 13. Jammer Matched Filtered Signal S 2 (t,u m ) (From [19]). 


A hard-limit (threshold) of -20 dB is applied onto S 2 (t,u m ) to chop off the peaks 
of the jammer signal, and the threshold output is denoted as S 3 (t,u m ). The hard-limited 
data is then transformed with a jammer inverse matched filter with the jammer pulse 
sequence so that the final signal, S 5 (t,u m ), contains the information only on the real 
target. The mathematical expressions are 

S 4 (t, u m ) = S 3 (t, u m ) * P m _, (t ) (4.6) 

S 5 (t,u m ) = S A (t,u m )*P;X-t) (4.7) 

An example picture of the resultant S 5 (t,u m ) is shown in Figure 14. The strength of the 
jammer signal is found to be approximately 40 dB less than that of the real target return. 
Although the algorithm can be further performed in azimuth domain for better 
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performance, as described in [20], the example in Figure 14 suggests DRFM penalizing 
algorithm in suppression of the DRFM false-target jamming. 
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Range Reconstruction Using Penalization & Radar Matched Filter 
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Figure 14. Jammer-Penalized Output Signal S 5 (t,u m ) (From [19]). 

C. CROSS-TRACK INTERFEROMETRY 

The dual-aperture antenna cross-track interferometer is a method, separate from 
those discussed above, that is proposed to locate the false target generated by the 
deceptive jammer. This method utilizes the phase infonnation of the image, pixel-by- 
pixel, and a spatial cancellation algorithm to separate false targets from real targets. 
Figure 15 shows a flow diagram of the dual-aperture cross-track interferometry scheme, 
in which it is easy to identify that the spatial cancellation is perfonned after the synthetic 
image is generated. Since this method does not involve complicated signal processing, it 
has great possibility to be applied on ISAR for the DIS false-target recognition, although 
it was originally proposed for SAR. 
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Antenna A 
(Receive-Only) 


Antenna B 
(Transceive) 



Figure 15. Flow Diagram of the Dual-Aperture Cross-Track Interferometry Scheme. 


To execute the cross-track interferometry method, two identical aperture antennas, 
SAR A and SAR B, must be aligned in parallel with a fixed displacement, d The 

geometry of the dual-aperture antenna is shown in Figure 16. 
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Figure 16. Geometry of the Dual-Aperture Antenna Cross-Track Interferometry 

(from [15]). 

SAR B operates in transceive mode, while SAR A is in a receive-only mode. 
Therefore, two synthetic images are obtained after the processing, namely image A and 
image B, respectively. Consider that a reflective surface, or a scatterer, is shown at 
location P on both images, and the slant range between SAR B and P is represented by r . 
There is a phase value with respect to the scatterer at P for both images, and the 
difference between those two phase values is [14] 

¥ t , = V, ~ V, = ^(AP - BP) * _ 2dx ) (4.8) 

Ar Ar 

where x p is the x-axis coordinate of P and A is the wavelength of the antenna radiation. 

Assume a false target is generated by a platform with a DRFM jammer located at 
P, and the x-axis coordinate of the false target is denoted as x p . Thus, according to (4.8), 

the phase difference at which the false target is located between two images can be 
expressed as 
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(4.9) 


■ i i > 

Yab = Wa -Yb~ “ 2 dx p ) 

Ar 

where r is the slant range between SAR B and the false target location. For simplicity, 
the false target is assumed to be in an extend such that r « r . Then, (4.9) can be rewritten 
as 

> i i 0 JT 

Yab = y'a -y'b -2 dx ) (4.10) 

Ar 

A weighting coefficient, w mn , is chosen to have the value of e ~ WAB , which is related to x p . 

This weighting coefficient is applied to the imaging output of SAR A, so that the 
nonnalized signal of each image pixel on image A can be expressed as 

S A = S A e i¥AB = e i¥A e j¥AB (4.11) 

while the nonnalized signal of the same image pixel on image B can be expressed as 

S B = e w " (4.12) 

The spatial cancellation block in Figure 15 essentially calculates the difference of 
the signal of the same location pixel between both images. The output of the spatial 
cancellation block is 


S M = S A e 


-JVab 


~S B =e 


JYa „~JYab 


-e 


JVb 


(4.13) 


Thus, if there is no false target and the image pixel, P , contains only the real target, the 
spatial cancelation output is 


c _ o p-iViB _ v 
^AB ~ ^A e ^B 

= Q iVA e~AV A -¥ B ) _ e W B 
= e J¥s _ £ Wb =0 


(4.14) 


If, on the other hand, the same image pixel contains both the real and false-target returns, 
the spatial cancelation becomes 
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(4.15) 


S AB =(S A +S A )e^-(S B +S B ) 

= ( e» A + e i¥A - (e i¥B + e j¥B ) 


= e 


Wb 


?K¥b-¥a+¥a-Vb) _ | 


= e 


J¥b 


..2nd ■ 


-1 


-0 


for the case x =x . Thus, from (4.14) and (4.15), it can be concluded that under two 

circumstances the output of the spatial cancelation on a pixel of an image is cancelled. 
Either there is no false target, or there is a false target and such false target lies on top of 
the real target within the same pixel of the synthetic image. In both cases, the image 
originally shown in the pixel is filtered. 

There is one condition left for discussion, i.e., a false target that lies at a certain 
extent from the jammer platform, or the real target. In this case, x p ^ x , and the output 

of the spatial cancellation is not 0. As a result, the signal of the real target is filtered while 
that of the false target remains on the image. After the completion of the process for 
every pixel on the image, the false target is portrayed and analyzed. The filtered real 
target is not completely lost; rather, it remains on the output image after the cancellation 
but at a weaker energy level, due to the inherent calculation errors in phase throughout 
the synthetic imaging processing. 

D. CHAPTER SUMMERY 

Different strategies have been considered to recognize, or even defeat, the false- 
target signal generated by the DIS on a DRFM platform. Three different methods based 
on the continuity of the target image, pulse diversity, and cross-track interferometry, 
respectively, are described in this chapter. The fundamental theory and mathematical 
expressions are also included within each individual discussion. 

Among these three methods, the methods based on the pulse diversity and cross¬ 
track interferometry are designed to reject the DRFM-based jammer either before the 
image compression processing or after the range-Doppler image is generated. In other 
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words, there is no better way for current DIS design to generate the false-target image on 
radars capable of these two EP strategies. 

Discriminating the DIS false target through determining the continuity of the 
target image is a method that takes advantage of the deficiency of the DIS design in the 
gap produced in cross-range domain. If a modification can be developed for the DIS 
design to compensate for such deficiency, which makes the DIS false-target image more 
continuous in cross-range domain, the false-target image will be more realistic and 
resistant to this discrimination method. 

In the next chapter, the addition of sea clutter to the DIS model is discussed in an 
effort to counter the image discontinuity EP technique used by ISAR. 
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V. THE DIS DESIGN WITH A SEA CLUTTER MODULE 


This chapter presents a modified DIS design incorporating a new sea clutter 
module. The purpose of adding the sea clutter module into the DIS design is to add a sea 
clutter effect to the false target, so that the interaction between the reflective surface of 
the target object and the background clutter is emulated, which makes the false-target 
image more realistic. It is shown that the Doppler shift derived from the sea clutter can be 
used to ease the discontinuity in the cross-range domain as described in previous chapters. 

The chapter starts with an introduction on theoretical methods that characterize 
the behavior of sea clutter, followed by an illustration on how the sea clutter module 
generates the sea clutter profiles using the DIS signal processing. 

The MATLAB code is composed of three major files and some minor files for 
parameter tracking purposes. Those main files must be executed in the order of 
extract_para_ship_v5, mathost_v5, and simhwchk_v5 . In extract_para_ship_v5, the 
modification coefficients for both DIS and the sea clutter module are obtained for later 
false-target generation. In mathost_v5, the DRFM samples used by DIS and the 
reference signal fed at the ISAR receiver are simulated. Finally, the DIS signal 
processing and ISAR image compression for producing a target range-Doppler image are 
simulated in simhwchk_v5. The simulation result of the modified DIS is provided and 
compared with that of the original DIS design. 

A. INTRODUCTOIN OF SEA CLUTTER MODELING 

Observations of sea clutter are normally associated with particular characteristics 
of the sea surface and the enviromnent, such sea waves, clutter reflectivity, and wind 
speed. Because of the dependency on various complicated environmental parameters, it is 
not an easy task to establish a model that precisely describes the behavior of sea clutter in 
detail. It is not always clear which environmental parameter is the most critical, and even 
if determined, it is often too difficult to measure such parameters accurately through the 
real sea conditions. 
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1 . 


Gaussian Clutter Model 


A complex Gaussian process is often developed to model the sea clutter in a low- 
resolution maritime radar system. When the radar illuminates a collection of scatterers (or 
clutter), each scatterer makes a contribution to the returned field. Each of these individual 
scatter fields is represented in a complex form that is known as the I and Q components. 
Consequently, each of the complex fields can be summed up to give a resultant vector to 
form the total scattered field, E , as follows 


N 

E = a „ and 

n =1 

(5.1) 

1 i<l 

a „ =d°n e 

(5.2) 


where a n represents the individual scatterer, N represents the number of scatterers, <J n 
is the scatterer cross section value, and (f) represents the scatterer phase. 

For a low-resolution radar system, the resolution cell is large enough to cover a 
large number of scatterers, which are considered as evenly distributed and decorrelated 
when illuminated by radar. Based on such assumptions, it is possible to develop a 
statistical Gaussian model, with a probability density function (pdf) for the clutter as 
follows [22,23] 

P(E I ,E 0 ) = — exp (-(£; + E\ ) / x); -^<E n E 0 <^ (5.3) 

nx 

The corresponding pdf of the envelope and intensity of the clutter signal can be expressed 
as 

2/7 I - 

P(E) = —exp (~E 2 /x); E = ME] + El)-, 0<£<oo (5.4) 

x v 

1 

P(z) = — exp(-z/x); z = E~; 0<z<oo (5.5) 

z 

where x is the mean intensity of the clutter return signal. The intensity can be directly 
linked to the received power in the radar equation in the following form 
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(5.6) 


P t G t G r A 2 cr° 

(4 7rfR 4 L 

where P t is the power of the signal transmitted by the radar, G t and G r are the gain of 
the radar, A is the radar wavelength, R is the distance between the radar and the jammer, 
<j° is the normalized clutter RCS, and L is the overall loss. 

Notice that, in (5.2), the phase of each scatterer is assumed to be constant, and 
thus, the Gaussian model in (5.3) carries a uniform distribution in phase. If the phase 
values are different, then this leads to a non-uniform distribution in phase for the 
Gaussian model, and is described in [24], 

2. Compound K-Distribution Model 

Although the Gaussian model is reasonable to represent the behavior of the sea 
clutter for calculating the perfonnance of low-resolution radar systems, the same model 
fails if applied to a high-resolution system, mainly due to the difference in resolution and 
the number of scatterers covered by a resolution cell. In a high-resolution radar, the sea 
surface has more resolution cells, such that some of the larger-scale structures of the sea 
surface can be further decomposed into several smaller structures. The decomposition 
makes the clutter much more variant than what is seen in low-resolution radar, i.e., the 
speckle-like and spike-like clutter. In this case, the assumption that supports the Gaussian 
model may not hold anymore, and a realistic non-Gaussian clutter model is needed. 

As mentioned earlier, the structure of the sea surface is very complex. Many 
independent small-scale structures within a single high-resolution cell can be recognized 
as even-distributed speckle-like clutters and described locally by the Gaussian clutter 
return, as shown in (5.5). However, those structures are possibly modulated by other 
more slow-changing large-scale structures at the same time, which implies a need to 
modify the local power derived from the Gaussian model. 

In order to take the dependency of the mean local power x into account, the pdf 
of the envelope of the clutter returns of the Gaussian model in (5.4) can be rewritten as 
[ 22 ] 
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0 < E < oo 


(5.7) 


P(E | x) = —exp(-fs 2 / x); 
x 

Here, the mean local power itself is a random variable with a pdf P(x) . By analyzing a 

sufficiently large quantity of data, it is found that the gamma distribution is best fit to the 
measured data, and thus 

TV 

P (x) =-X 1 '" 1 exp(-bx); 0 < x < oo (5.8) 

T(v) 

where the gamma distribution is described by a gamma function, T(»), a scale parameter, 
b , and a shape parameter, v. 

T The pdf of the envelope of the non-Gaussian clutter is obtained by integrating 
the local Gaussian clutter returns with mean local power x over P c (x) [22,25] 


/•oo 

P(E) = j P(E | x)P(x)dx 

/•CO | 

= 2iiJ —exp(E\x)P c (x)dx 

X 


(5.9) 


When P(x) is replaced by the result from (5.8), the pdf in (5.9) can be rewritten as [22] 


P(E) = f x 1 ' 2 exp(-hx - E 2 / x)dx 

r(v') Jo 


4E v b (v+l)l2 

r(v) 


K,_,(2Esfb) 


and the corresponding pdf of the intensity becomes 


2 / 1 ( v+1 ) J y - l ) /2 

P{z ) =-2_-(2 yfbz); where z = E 2 

T(v) 


(5.10) 


(5.11) 


where K function represents the modified Bessel function. From (5.7) to (5.9), a general 
model with compound pdf is defined to emulate the non-Gaussian clutter (sometimes 
called the K-distribution model due to the utilization of the K function). 
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he K-distribution model, as described, was developed by K.D. Ward to simulate 
sea clutters when observed by a high-resolution radar. This model suggests that the 
amplitude of the K-distributed clutter within a resolution cell can be generated by a 
combination of two statistical random numbers, the Gaussian random number and the 
gamma random number [22,23]. 

Except for the speckle-like clutter, there are events where the discrete spike-like 
clutter (i.e., the burst) occurs within a radar resolution cell. The K-distributed model is 
further modified to incorporate the spike, which results in what is known as the KA 
model. The corresponding pdfs for the Gaussian and gamma distributed components 
become [22] 


CO 1 

P(z 1 *) = y exp 

f -z ' 

P, (m); 0 < z < oo 

(5.12) 

m=0 x + ma sp 

i 

[x + m( 7 sp j 



P(x) = , . x 
cV ’ r(.) 

1 exp(-bx); 

0 < x < oo 

(5.13) 


where a sp represents the local mean spike power, m represents the number of spikes in a 
resolution cell, x represents the local mean speckle power, and P m (m ) represents the 

probability of m spikes occurring in a resolution cell based on the Poisson distribution. 
In this model, it is assumed that the mean number of spike-like clutter in each resolution 
cell is N , and TV <sc 1, so that P m (/??) is defined as [22] 

^.(0) = 1-JV 

^(1) = V 

p, (»< s 2) = o 

which, in other words, the probability for more than two spikes occurring within a range 
cell is zero. 
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3. Doppler Spectrum Based on K-Distributed Model 

The structure of the sea surface can sometimes be complex due to the inconsistent 
variation of the compound component of the sea structure over an observation time 
interval T . Some structures may vary in less than a second while others on the order of 
many seconds, and different choices of T may suggest different results when analyzing 
the sea clutter. In practice, such T is carefully chosen to be at least as large as the time 
taken for the cycle of slow-varying large-scale structures that contain the modulation of 
some small-scale structures in the local speckle power. 

In the K-distribution model, the gamma-distributed local power x is randomly 
determined over a certain time interval in which x has not changed perceptibly. Thus, 
such v can be considered as the integral of the power spectrum over all frequencies 
within that interval. As a result, the K-distributed model can be transformed into a model 
in frequency domain that describes the fluctuating Doppler spectrum as the product of the 
local power x and a suitable function of unit area [22] 

S(«|x) = xS(«) (5.14) 

oo 

J S(co)dco = 1 (5.15) 


where S ( co ) itself can be modelled as a random Gaussian process similar to (5.4) in the 
following form 


5 (®) 


exp (co - co 0 ) 2 jla 2 j 
4 2kg 1 


(5.16) 


From (5.15) and (5.16), it is easy to understand that this model is based on the 
simple assumption that the frequency shift in all Doppler bins is identical. However, 
according to the observation on sea structures, larger frequency shifts usually occur on 
those of spike-like clutters. In order to make the power spectrum model more applicable, 
other than the speckle-like clutters, it is further associated with the spike-like clutters. 
Thus, the Doppler spectrum consists of two Gaussian components as below 
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(5.17) 


c / \ _ x 01 S 0 (®,^o 5 °) + ^^i(^ cr P Q i) 

o ^ CO J — Xj- 


x 0 + x l 


S ( co , cr,Q) 


exp(-(<y-Q) /2cr 2 ) 
\]l no 2 


(5.18) 


where x is the gamma-distributed local power that modulates the averaged clutter power 
over the frequency spectrum and Q is the frequency shift from the center. 

In (5.17), the first Gaussian component with a zero frequency displacement and 
relatively small power distribution, S 0 (<z>,cr 0 ,0), represents the speckle-like clutter from 
small-scale structures, while the second one with a frequency displacement from the 
center and a relative significant power distribution, represents the spike¬ 

like clutter from large-scale structures. 

B. CLUTTER GENERATION PROCESSING 

The sea clutter module developed emulates the sea clutter in terms of the 
amplitude and frequency displacement in the form of phase representation. Both the 
amplitude and phase information are then modulated into the gain and phase coefficients. 
Thus, new gain and phase coefficients are produced by changing (3.9) and (3.13) into 

</> T (r,n) = angle [real j T (r,n) + Clutter (r,n)],imag{r (r,n) +Clutter (r,n) }) (5.19) 


7 

v (r,n) + 

Clutter (r , n) 


max 

( T(r,n ) 

+ 

Clutter (r, n) j 


(5.20) 


Figure 17 shows the block diagram of the modified DIS design integrated with the 
sea clutter module. 
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Figure 17. DIS Block Diagram with Sea Clutter Module. 

1. Gain Coefficient Generation 

The generation of the gain coefficient of the sea clutter module is similar to what 
has been described in Chapter III, with the target-profde amplitude replaced by the 
amplitude of the sea clutter, which is generated as follows. 

According to the K-distributed model, the sea clutter can be statistically expressed 
by a combination of the gamma and Gaussian components with an additional clutter 
modulation in power on the gamma component [23]. A block diagram depicting the 
processing of generating the K-distributed clutter is shown in Figure 18. 
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Figure 18. Block Diagram of K-Distributed Clutter Generation Process. 


The RCS value for each individual clutter is randomly assigned based on an 
empirical RCS model developed by the Georgia Institute of Technology (GIT). This 
model is also known as the GIT model. This GIT model determines the relationship of 
the normalized clutter RCS versus grazing angle and sea state through wave ta nk 
experiments for radars of frequency ranging from 1 to 100 GHz. From the trend observed 
in the experimental data sets, the clutter RCS for horizontal and vertical can be 
approximately expressed as [25] 
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65.91 + 101084 W‘UA4.) 


(5.21) 


0 

= 


<jy = cr^ — 1.05 In (/z ave + 0.015) +1.09 In (/l) +1.27 In + O.OOOl) + 9.7 (5.22) 

where 1 is the radar wavelength, ^ is the grazing angle, 4 is the sea state dependency 
parameter, A u is the wind direction parameter, and A i is the multipath interference 
parameter. 

Figure 19 shows the experimental result at 9.5 GFIz for sea state 1 through 5, in 
which most of the RCS values fall into the range between -30 dB and -40 dB for both 
polarization at low grazing angle (less than 10 degrees) [22], Thus, in the simulation, the 
clutter RSC value is randomly chosen in the range between -20 to -40 dB to be more 
flexible to various sea state conditions and antenna polarizations. 




a. Vertical polarization b. Florizontal polarization 

Figure 19. Normalized RCS Versus Grazing Angle from GIT Model for Sea States 1 

to 5 (in the Order of Bottom to Top) (From [22]). 


According to (5.1) and (5.2), the number of the clutter is directly related to the 
intensity of the gamma-distributed clutter. Thus, in the sea clutter module, the number of 
the clutter within each range-Doppler cell is assigned by a random integer, which is 
intended to reflect the environmental variation. The mean local power, x , is then 
obtained through a MATLAB built-in gamma random function with the consideration of 


50 








the clutter RCS and distance between the false target and the radar. By the same token, 
the Gaussian-distributed clutter can also be obtained through a MATLAB built-in 
Gaussian random function based on x , and thus the final K-distributed clutter is 
generated. 

According to (3.7) and (3.8), the amplitudes of this K-distributed clutter for each 
range-Doppler and range cells are obtained as follows 

Clutter(r, d, n) = A clutter (r, d)e^ 1 2 * / * w(r,<0 " ra/ (5.23) 

Nd{r) 

Clutter (r,n) = ^ Clutter(r,d,ri) (5-24) 

d =1 

where A clutter (r,d) is the amplitude of the K-distributed clutter and f dutter (r,d) is the 

corresponding frequency displacement discussed in the next section. The gain coefficient 
can be obtained through the normalization process through (5.20), (5.23), (5.24), and 
Table 1. Figure 20 shows the normalized amplitude range profile for the clutter generated. 
A plot of the function Y - MR 2 marked in red is also included for a comparison purpose. 
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Figure 20. Normalized Amplitude Range Profile of the Sea Clutter. 
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With a closer examination on both Figure 20, it is clear that the amplitude of the 
clutters decreases dramatically in a behavior similar to that of the function Y -V R 2 . 
Such similarity is corresponding to the relationship between the power of the clutter and 
the distance, x oc 1 / R 4 , as described in (5.6). It implies that the amplitude of the clutters 
generated is appropriate for the production of the gain coefficient for the clutter. 

2. Phase Coefficient Generation 

The generation of phase coefficient in the sea clutter module is based on the same 
schemes to generate the DIS phase coefficient. The first step is to assign a frequency 
displacement for each clutter contained within a range-Doppler cell, so that the phase 
information of the clutter is obtained by taking the angle function on the clutter signal. In 
order to cooperate with the DIS and DRFM system, such phase information must also be 
quantized to fit the DRFM architecture. 

Although a theoretical model for a power spectrum is presented in the previous 
section, the determination of the frequency displacement for the sea clutter through this 
model is far complicated. An empirical model based on the analysis of the data collected 
from radar surveillance operations and wave tank experiments is adopted instead to 
simplify the processing of random frequency displacement generation instead. 

According to the experimental data from [22], [23], and [26], the empirical model 
is suggested to estimate the frequency displacement of the sea clutter for low grazing 
angles. It is found that the Doppler velocity at peak of the upwind spectrum appears to be 
approximately determined by the velocity of the wind speed for both polarizations in the 
following linear forms [26] 


Vd verbal ~ 0.25 + 0.1817 or 

(5.25) 

V D horizontal ~ 0-25 + 0.2U 

(5.26) 


where U represents the wind speed in meters per second. Since there is little difference 
between both polarizations, (5.26) is used to calculate the Doppler velocity value, V D , in 
the sea clutter module, and the frequency displacement can easily be determined by 
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(5.27) 


2V 

f _ D 

J D * 

radar 

here A radar represents the radar wavelength depending on the radar frequency. Table 2 
shows the approximate upwind frequency displacements measured by radars of 1 and 10 


GHz for sea state 1 to 4 by using (5.26) and (5.27). 



Wind Speed 

Wind Speed 

V D 

fD 

fD 


(kn) 

(m/s) 

-0.25+0.2U 

a=0.3m) 

a=0.03m) 

Sea State 1 

0-6 

0-3.0867 

0.25-0.5873 

1.66-5.78 

16.6-57.8 

Sea State2 

6-12 

3.0867-6.1733 

0.8673-1.4847 

5.78-9.898 

57.8-98.98 

Sea State3 

12-15 

6.1733-7.7166 

1.4847-1.7933 

9.898-11.955 

98.98-109.55 

Sea State4 

15-20 

7.7166-10.2888 

1.7933-2.0378 

11.955-15.385 

109.55-153.85 


Table 2. Upwind Frequency Displacement under Various Sea States 


When the wind turns from upwind to crosswind and to downwind, the frequency 
displacement shifts correspondingly from more positive displacement to zero and to 
negative value. Thus, the frequency for a clutter is randomly assigned as either positive or 
negative to represent the effect under the wind direction. An example of the code 
segment used for generating the frequency displacement for sea state 4 is shown below to 
illustrate the implementation of this frequency empirical model. 

%GENERATION OF FREQUENCY DISPLACEMENT 

if rand(l)<=0.5 %assign random wind direction 

clutter_freq=randi([-l 10,-100],[l,l])+2*randn(l,l); %negative frequency for downwind 
else 

clutter_freq=randi([ 100,110],[l,l])+2*randn(l,l); %positive frequency for upwind 


From the code segment listed above, the Doppler frequency of the clutter is 
chosen from a frequency spectrum constructed by a statistical Gaussian random function. 
The mean of the Gaussian random function is determined by the sea state, the radar 
wavelength, and the wind direction. The Gaussian standard deviation, on the other hand, 
is determined by the choosing a value of 2 in order to include 95% of the frequency 
within the spectrum. It is also assumed that all clutters have a non-zero Doppler for a 
more general module. 
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Clutter observed by a 10GHz radar under sea state 4 is generated and a pulse-to- 
pulse correlation is performed as shown in Figure 21 and demonstrates the correct drop 
off rate expected. The Doppler amplitude spectrum for range cell 1,5, 10, and 20 are 
examined as shown in Figure 22 to Figure 25. The amplitude frequency spectrum is 
obtained by taking the Fourier transform on the clutter signals contained within one range 
cell. Therefore, the characteristic of decreasing in amplitude is reflected in the figures 
when they are compared in the amplitude domain. 



Figure 21. Pulse-to-Pulse Correlation on the Clutter Generated. 
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Figure 22. Amplitude Frequency Spectrum at Range Cell 1. 
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Figure 23. Amplitude Frequency Spectrum at Range Cell 5. 
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Figure 24. Amplitude Frequency Spectrum at Range Cell 10. 



Figure 25. Amplitude Frequency Spectrum at Range Cell 20. 
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According to Table 2, the mean of the Gaussian random function is chosen 
approximately either between 100 and 110 or between -110 and -100. As a result, the 
center of the corresponding frequency spectrum is most likely between 50 and 100 in 
both positive and negative. 

The spike-like clutters with stronger energy level can be recognized in Figure 22, 
and their frequency values are found to be slightly different. Such difference could 
possibly be due to the correlation with the speckle-like clutters. The speckle-like clutters, 
on the other hand, have the frequency values spread through the spectrum with weaker 
energy. This suggests that such clutters travel for a certain amount of time and do not 
change if observed within a relatively short period of the time. In either case, the 
frequency displacement generated for the clutter with this empirical model is appropriate 
to represent the Doppler shift of a random clutter. 

C. SIMULATION RESULT 

1. Evaluation on Performance of the Modified DIS Design 

Figure 26 shows two ISAR images from the simulations based on an identical 
false-target profile of a ship with 32 by 256 range-Doppler cells. The image for the 
profile generated by the original DIS refers to Figure 26a, while that for the profile 
generated by the DIS incorporated with the sea clutter module is in Figure 26b. Although 
both images are discontinuous in the Doppler domain, which is inevitable due to the 
quantization of the phase into integer representation, it is clear to see that Figure 26b 
shows better continuity in both range and cross-range domains than Figure 26a. 
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a) Original DIS b) DIS Modified with Sea Clutter Module 

Figure 26. ISAR Image Simulation for DIS and Modified DIS. 

Also, notice that in the Doppler domain of the image in Figure 27b, there are 
some background noises forming in a shape similar to the ship around the structures with 
larger frequency shifts like the mast. This is the effect due to the correlation between the 
false-target and the sea clutter, which is the purpose behind incorporating the original 
DIS design with a sea clutter module. 

2. The Impact of Radar Wavelength 

According to (5.27), the frequency displacement generated by the empirical 
model has an inverse relationship with the radar wavelength. The purpose of this 
simulation is to examine how the radar wavelength could have impact on the empirical 
model and the modified DIS as well. 

Figure 27 shows two ISAR images from the simulations based on the false-target 
profile of a ship that rotates in negative frequency shifts against the radar with 32 by 256 
range-Doppler cells. These negative shifts make the target image upside down. The false- 
target profiles used for generating these two images are identical except for the frequency 
displacement generated by the sea clutter module. The frequency displacement chosen for 
Figure 27a is based on a radar wavelength of 0.3 meters for sea state 4, which, according 
to Table 2, has either a positive value approximately between 12 and 16, or a negative 
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value approximately between -16 and -12. The frequency displacement chosen for Figure 
28b, on the other hand, is based on a radar wavelength of 0.03 meters and has either a 
positive value approximately between 120 and 160, or a negative value approximately 


between -160 and -120. 
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Figure 27. ISAR Image for Modified DIS at Different Radar Wavelength. 


With a closer inspection on both images, other than some variances of the 
background noise distribution, the majority of the false-target images are in great 
similarity. In other words, the frequency empirical model used to generate the frequency 
displacement has little relationship with the radar wavelength (or frequency) on which 
ISAR usually operates. 

3. The Impact of Sea State 

The sea state is one of the most critical environmental variables that can greatly 
affect the RCS characteristic of the sea clutter. As suggested from Table 2, the higher the 
sea state, the more complicated and rough the sea surface. In the models used for the sea 
clutter module, the generation of the clutter RCS and frequency displacement for the 
clutter is also involved with the sea state, and thus, this simulation is performed to 
examine if the false-target profile with the randomly assigned clutter RCS and frequency 
displacement produces a reasonable image in the same manner as varied sea state. 
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Figure 28 shows ISAR images from the simulations based on a false-target profile 
of a ship that rotates in negative frequency shifts against the radar with 32 by 256 range- 
Doppler cells for sea state 1 through 4. These negative shifts make the target image 
upside down. The intensity of the ship body and clutter is identical for all four images. 



a) Sea State 1 b) Sea State 2 




5 10 15 20 25 30 35 ■ 5 10 15 20 25 30 35 

_Down Range Cells_ U _Down Range Cells_ 


c) Sea State 3 d) Sea State 4 

Figure 28. ISAR Image for Modified DIS at Different Sea State. 

Such results suggest there is no tremendous difference in the clutter RCS with 
respect to various sea states for a low grazing angle, which corresponds with the 
observation from the experimental data, as shown in Figure 19. As far as the frequency 
displacement is concerned, the amount of clutter with a large frequency displacement 
increases as the sea state gets higher. Such clutter can be referred to as spike-like clutter 
caused by the increasing wind speed at higher sea state. Hence, the empirical models used 
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for generating the clutter RCS and frequency displacement are able to reflect a realistic 
behavior of clutter in accordance with the change in sea state. 

D. CHAPTER SUMMERY 

In the beginning of this chapter, different theoretical models used to statistically 
characterize the behavior of the sea clutter in amplitude and frequency displacement with 
respect to critical environmental parameters are presented. Next, an illustration is 
presented on how the sea clutter module generates sea clutter profiles compatible with the 
DIS signal processing to enhance the false-target signal through empirical models based 
on the theoretical methods. Finally, a series of discussions on the simulation results are 
provided to verily the perfonnance of the novel modified DIS design and the 
reasonableness of the methodology for random clutter generation performed by the sea 
clutter module. The next chapter provides the thesis conclusions and concluding remarks. 
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VI. CONCLUSION AND REMARKS 


This thesis first introduces the ISAR imaging process followed by a discussion of 
the theoretical operation of electronic attack against ISAR, as well as the electronic 
protection of ISAR from the DRFM-based DIS false targets, through research. Next, a 
modified DIS design that incorporates a sea clutter module is introduced as an improved 
EA technique against ISAR. In order to generate reasonable sea clutter, several statistical 
sea clutter models are described, as well as the random processes used to generate the 
amplitude and frequency displacement coefficients that represent the random sea clutter. 
Finally, the simulation of the modified DIS design is executed to evaluate the 
performance of the new design. Several observations and remarks are concluded 
throughout the thesis and are described as follows. 

A. CONCLUSION 

The algorithm used to discriminate DRFM-based false targets with the pulse 
diversity technique could suppress the DRFM-based jamming signal up to -20 dB lower 
than the true radar returns and thus could easily reject the jamming signals sent from 
DRFM-based jammers. The DRFM penalizing algorithm, on the other hand, could 
suppress the jamming signal up to -40 dB lower than the true radar returns on the premise 
that the characteristic of the jammer is known beforehand. However, there is a great 
difficulty to obtain such information. Thus, the algorithm with the pulse diversity 
technique is much more flexible than the DRFM penalizing one even if it does not have 
the best performance against the DRFM-based jamming technique that covers the one 
derived from the DIS design. 

With two identical radars operating at the same time, one in transceive mode and 
the other in receive-only mode, the cross-track interferometry scheme is able to locate the 
false target within a certain distance from the DRFM jammer platfonn through a spatial 
cancellation process on the output image of the two radars. Although both radars receive 
the false-target signal, it is the displacement between them that produces a delay, or a 
phase difference, that allows the cancellation process to partially locate the false target. In 
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order for the spatial cancellation to function properly, the false image must be carefully 
generated so that it has the slant range with respect to the transceive radar approximately 
equal to that of the jammer platform. Such limitation reduces its flexibility against the 
DRFM-based false target. Nevertheless, the cross-track interferometry algorithm can be, 
in general, considered as an effective method to counter the DRFM-based false-target 
images, including those generated by the DIS. 

Determining the continuity of the target image in cross-range domain is one of the 
most convenient methods to discriminate the DIS false target since it takes advantage of 
the flaw of the DIS design in the possible gap produced during the digitization process on 
the phase information, including the phase sampling of the intercepted ISAR signal and 
the phase coefficient generated by the DIS microprocessor. Such errors lead to the 
discontinuity of the image in cross-range domain, which makes the images less likely to 
be a real target. The design of the DIS false-target profile also contributes to this 
discontinuity. The modified DIS design proposed in this thesis adds a sea clutter effect to 
the DIS false target so that the interaction between the reflective surface of the target 
object and the background clutter is emulated. In addition, a series of simulations prove 
that the random Doppler shift of the sea clutter derived from the modified DIS design can 
be used to ease the discontinuity in cross-range domain. 

B. RECOMMENDATIONS 

There is still room for the MATLAB software to be improved by replacing some 
of the current iterative loop structures into matrix operation so that the performance of the 
parallel operations in DIS is increased and the simulation time is reduced, as well. 

The amplitude and frequency displacement of the sea clutter is randomly 
generated through a Gaussian and gamma statistics model based on unified parameters, 
which may not be optimal. More research can be performed to identify a better setting for 
the parameters of the related statistic model. Hence, the amplitude and frequency 
displacement of the clutter can be assigned with an even more realistic value. 
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